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We have built a code to numerically solve the Riemann problem in ideal magnetohydro- 
dynamics (MHD) for an arbitrary initial condition to investigate a variety of solutions 
more thoroughly. The code can handle not only regular solutions, in which no interme- 
diate shocks are involved, but also all types of non-regular solutions if any. As a first 
applic ation, we explored the neighborhood of the initial condition that was first picked 
up by iBrio fc Wu ( 1988 ) and has been frequently employed in the literature as a stan- 
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dard problem to validate numerical codes. Contrary to the conventional wisdom that 
there will always be a regular solution, we found an initial condition, for which there 
is no regular solution but a non-regular one. The latter solution has only regular solu- 
tions in its neighborhood and actually sits on the boundary of regular solutions. This 
implies that the regular solutions are not sufficient to solve the ideal MHD Riemann 
problem and suggests that at least some types of non-regular solutions are physical. We 
also demonstrate that the non-regular solutions are not unique. In fact, we found for the 
Brio fc Wu initial condition that there are uncountably many non-regular solutions. This 
poses an intriguing question: why a particular non-regular solution is always obtained 
in numerical simulations? This has important ramifications to the discussion of which 
intermediate shocks are really admissible. 

PACS codes: 



1. Introduction 

The Riemann problem for a system of Ist-order, hyperbolic partial difi^erential equa- 
tions is an initial value problem, in which two distinct constant states are separated by 
a discontinuity initially. By hyperbolic we mean that all eigenvalues of the Jacobian ma- 
trix that characterizes small perturbations are real. It is well known that the Riemann 
problem in ideal MHD has in general non-unique solutions for a given initial condition 
even if one discards manifestly unphysical ones such as those including rarefaction shock 
waves, in which entropy is decreased. This is in sharp contrast to the hydrodynamical 
counterparts and is due to the facts that the system of equations that describes ideal 
MHD is not strictly hyperbolic, i.e. some of the eigenvalues of the Jacobian matrix are 
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coincident with each other, and that the c haracteristic fields are n either hnear nor gen- 
uinely non-linear (see §[2]and e.g. Laxlll957 : Jeffrey fc Taniuti 19641 : Polovin fc Demutskiil 
19901) . In fact, the MHD Riemann problem has a far greater variety of solutions including 
so-called intermediate shock waves. Although these intermediate shocks are deemed to 
be unphysical in some text books , some interplanetary experiments have reported the 
detections of intermediate shocks ( Chaol[l995t iFeng fc Wang|[2Q08l: iFeng et aLll2009l ) and 
they are actually realized in almost all numerical solutions. Their reality is hence still an 
issue of controversy as will be described below more in detail. 

Solutions of the Riemann problem consist of discontinuities and centered self-similar 
simple waves in general (see §[2] for more details). Various shock waves appearing in the 
solutions of Riemann problem are usually characterized by the flow speeds upstream and 
downstream in the shock-rest frame: in ideal MHD, for example, the fast shock has a 
flow velocity that exceeds the fast-wave velocity (termed super-fast hereafter) upstream 
whereas the downstream flow velocity is sub-fast, i.e. lower than the fast-wave velocity, 
but super- Alfvenic, i.e. higher than the Alfven-wave velocity (see §[3]for the definitions of 
these wave velocities); assigning the numbers 1, 2, 3 and 4 to the super-fast, sub-fast and 
super- Alfvenic, sub- Alfvenic and super-slow, and sub-slow, respectively, we designate the 
fast shock as 1 — > 2 shock; other shocks are referred to just in the same way. 

The intermediate shock waves in ideal MHD have coplanar configurations, i.e. the 
upstream transverse magnetic field is anti-parallel to the downstream one and they are 
characterized by the flow velocities that are super-Alfvenic upstream and sub- Alfvenic 
downstream. They are normally assigned to one of the following four types: 1 — ?> 3, 1 — > 4, 
2 — > 3 or 2 — > 4 shock. There are some boundary types as well: for example, the 2 — > 3, 4 
shock is a boundary type that has a downstream flow speed that is equal to the slow 
speed as indicated by its name. As mentioned above, these intermediate shock waves 
are considered to be unphysical i n some textbooks (see e.g. iLax 1957; Je ffrev fc Taniuti 
19641 : Polovin fc Demutskii 1990[ ) because they do not satisfy the so-called evolutionary 



conditions, which require the existence of perturbed states of the same type; for ordinary 
shock waves (such as the fast shock in ideal MHD), one of the characteristics runs into the 
discontinuity from both sides; this is not the case for some of the intermediate shocks; for 
instance, the 1 — >■ 3 intermediate shock is a shock wave whose flow speeds in the shock- 
rest frame are super-fast upstream and sub- Alfvenic and super-slow downstream; hence 
the fast and Alfven characteristics go into the shock both upstream and downstream; in 
the case of the 2 — 3 shock, only the Alfven characteristic converges to the shock wave; 
however, another evolutionary condition on the linear independence of eigenfunctions is 
violated and it is classified as unphysical shock. 

In spite of the evolutionary conditions, the intermediate shocks are commonly observed 
as a stable constituent wave in num e rical solutions of the MHD Riemann problems 



19871 : iBrio fc Wij [19881 : Iwi] ll988&iy. Il99nl : |Wu fc Kenne]ll992h . One of the best known 
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cases is the one discussed bv iBrio fc Wul ( 19881 ). in which the numerical solution con- 
tains a so-called compound wave that consists of a 2 — >■ 3, 4 shock followed by a slow 
rarefaction wave. lWul (1987) also showed numerically that the 2 — > 4 shock- like structure 
could be obtained by steepening of a continuous wave, an indication that the resultant 
structure is physical. In addition he demonstrated that a perturbed states obta i ned for 
a slightly different initial condition retains a similar structure. Moreover. I Wul (19886) 
computed interactions between the 2 — )■ 4 shock-like structure and an Alfven wave and 
found that the structure remains for a while before it breaks to other waves. Thes e results 
indicate that the intermediate shock is stable to the perturbation in a sense. I Wul (|l9886r ) 
found that the shock-like structures observed in the computations with perturbations are 
time-dependent and do not satisfy the Rankine-Hugoniot condition in fact. He claimed 
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that these new types of shocks are the neighboring states for the intermediate shocks. 
In yet another simulation Wu ((1988al showed that rotational discontinuities, which are 
not intermediate, break into some other waves that include an intermediate shock. This 
is a suggestion that it is not intermediate shocks but particular type s of evolu tionary 
discontinuities that are unstable and unphysical. On the other hand, Wu (|l99Ch ) stud- 
ied the structures inside shock waves by introducing dissipations and found that they 
are not uniquely determined by the asymptotic states and there remain extra degrees 
of freedom. He concluded that the stability of the intermediate shocks to the interac- 
tions with Alfven waves is originated from these degrees of freedom. Incidentally, it was 
also demonstrated that all types of intermediate shocks can be produced by steepening 
processes in dissipative MHD. 

Complete ly opposite and equally convincing arguments were made bv lFalle fc Komissarovl 
(|l997ll200lh . They pointed out that if perfect planarity or coplanarity is imposed initially, 
the solution retains the symmetry and magnetic fields cannot rotate in the evolutions; 
then the problem is equivalent to solving the reduced system of MHD equations, in 
which it is assumed that magnetic fields are confined in a plane and, as a consequence, 
the Alfven waves do not exist; this changes the number of characteristics of the system 
and modifies the evolutionary conditions; in fact, the 180°-rotational discontinuity be- 
comes non-evolutionary whereas 1 — > 3 and 2 — >■ 4 shocks become evolutionary for this 
reduced system. They also demonstrated numerically that the intermediate shock does 
not emerge if coplanarity of the solution is broken by insertin g a thin layer within which 
the magnetic field rotates continuously. Barmin et al\ ( 1996[ ). on the other hand, found 
numerically that the compound wave breaks into a r otational discontinuity and a slow 
shoc k if the exa ct coplanarity is perturbed. Although iFalle fc Komissarov ( 2001 ') agreed 
within (|l990t ) that the temporary survival of 1 — )■ 3, 2 — )■ 4 and 1 — )■ 4 shocks in their 
interactions with Alfven waves is due to the non-unique shock structures, it was claimed 
that they should be regarded as transients. 

The evolutionary conditions themselves have b een reconsidered in the mean time (jHada 
19941 : Markovski] 1998 : Inoue fc Inutsuka 2007 ). These authors took dissipations into 



account, considering that the waves that emerge in dissipative MHD should be physical. 
With dissiptations shocks are no longer discontinuities and the analysis is facilitated. 
They found new modes that do not exist in ideal MHD and argued that on account of 
these modes the intermediate shocks are indeed evolutionary. 

Putting aside the reality of the intermediate shocks, iTorrilhonI (|2002l . 120031 ) studied 
when the intermediate shocks appear i n solutions of the MHD Riemann problem more 
in detail. In his paper (jTorrilhonl 120031 ) . he defined at first the regular and non-regular 
shocks: the former has A'^-t-1 incoming characteristics, where N is the number of equations, 
whereas the latter includes all the rest. The definition is a bit vague for those shocks that 
have a characteristic speed equal to the shock velocity, since it is not obvious in this case 
that the characte ristic is incoming or out going and, in fact, there has been a controversy 
in the literature ( Jeffrev fc Taniuti 1964J) . Although there is no explicit mention on this 
issue, the author of (iTorrilhonl l I2OO21 120031) regarded th em as non-regular We will follow 



this practice hereafter in this papei|jj. ITorrilhonI (|2003f ) referred to the solutions without 
non-regular shocks as regular solutions and all the others are called non-regular solutions. 
Believing that there is always a regular solution, he concluded that if a certain condition 



f More precisely, we define the regular shocks as the shocks that satisfy the evolutionary 
condition (see i )2.3[) . We treat those characteristics that have the same velocity as the shock 
velocity as incoming. All the other shocks are defined to be non-regular. The solutions are 
referred to as regular if they include no non-regular shocks and called non-regular otherwise. 
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is satisfied, there is a unique regular solution to the ideal MHD Riemann problem. Note, 
however, his analysis was not complete, since he assumed that only a left-going or right- 
going 1 — 3 or 2 — 4 intermediate shock exists and he ignored the cases in which the 
transverse magnetic field vanishes. Furthermore, he did not prove that there is always a 
unique regular solution. 

We address these issues in this paper, extending the approach in lTorrilhonI (j2002ll2003r ). 
We believe that the detailed study of particular solutions of the MHD Riemann problem 
together with their neighboring ones are also helpful to get some more insights into the 
reality of intermediate shocks. For this purpose we coded a program to produce both 
regular and non-regula r solutions for arbitrar y initial conditions. We have extended the 
algorithm proposed by iTorrilhonI (I2002L l2003l) so that we could handle all combinations 



of con stituent waves, including all types of intermediate shocks. Note that lAndreev et al 



( 2008[ ) are providing a code to solve the MHD Riemann problem that is available online. 
However, it can handle only the regular solutions. It should be emphasized that numerical 
codes to compute temporal evolutions are useless for the current purpose, since they 
produce a single solution for a given initial condition and we do not know how numerical 
dissipations that are inevitably inherent to numerical simulations affect the outcomes. 
Hence we need to construct solutions in all possible ways to satisfy the given jump 
condition. 

The parameter space of the MHD Riemann problem is vast (15 dimensions indeed) and 
exploring the entire space is almost impossible. We hence need to pay attention to those 
initial conditions and their neighborhoods, which appear of particular interest to us in 
the context explained above. In th is paper we i nvesti gated two cases: the first one is the 
initial condition first picked up bv lBrio fc Wu ( 1988[) that includes the compound wave 
in one of the solutions as mentioned earlier; we found that there is a regular solution in 
this case and that there are in fact uncountably many non-regular solutions as well; this 
poses an interesting question why the non-regular solution with the compound wave is 
always realized in numerical simulations; in the second case, we found that there is no 
regular solution; this provides a counter-example to the conventional wisdom that there is 
always a regular solution. Note that the regular solutions normally consist of fast waves, 
slow waves, rotational discontinuities and a contact surface with an exception that there 
exit so-called switch-off rarefaction waves that nullifies the transverse magnetic fields. 
Then rotational discontinuity does not appear in the solutions. 

The paper is organized as follows. In § [2l we briefly review the general theory of 
Riemann problem. In § [3l the shock waves and simple waves in ideal MHD, which are 
the constituent waves in the solutions of MHD Riemann problem, are presented. In § [H 
we describe how to construct various solutions that satisfy the jump conditions. In § [5l 
we apply the code to two particular initial conditions mentioned above and show both 
the regular and non-regular solutions. We discuss possible implications of the results in 
§|6l 



2. Theory of Riemann problem 

2.1. Hyperbolic equations and weak solutions 

In this section we consider general conservation equations in one spatial dimension. Since 
the Riemann problem assumes plane-symmetry, this is sufficient. Then the coupled equa- 
tions are expressed as 

du df ^ 
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where it is a vector of conserved variables and f{u) is a flux vector. The system (|2.ip is 
called hyperbolic if the Jacobian matrix, which is defined as 

dfi 

has N real eigenvalues \k{u) (k = 1,2,..., N) and corre sponding N linealy independent 
right eigenvectors, rk{u) fe.g. iFalle fc Komissarov 2001 ). In particular, if all eigenvalues 



are different from each other for any u, the system is called strictly hyperbolic. 

The /c-th characteristic field is called genuinely non-linear if the following condition is 
satisfied for any u: 

Vk ■ V^Afe ^ 0. (2.3) 

Here Vu denotes the operator, *(9/9ui, . . . ,d/duN), where *(. . .) stands for transposi- 
tion. On the other hand, The fc-th characteristic field is called linear if 

Tk ■ V^Afc = (2.4) 

for all u. Note that the characteristic field may be neither linear nor genuinely non-linear. 

In the system described by Eq. (|2.ip . discontinuities are ubiquitous. Indeed they may 
be spontaneously produced from continuous initial conditions by wave-steepening. Since 
Eq. (|2.ip is differential equations, these discontinuities cannot be treated as they are and 
some extensions are needed. The weak solutions of Eq. (|2.1I) are defined as the functions 
that satisfy the following equations instead of Eq. (12. ip in a finite domain D{x,t) (e.g. 
[jefl^rev fc Taniuti 1964): 

J w{x,0)-ct,{x)dx + At j da;f — .«+ — . /j =0, (2.5) 

where (f){x) is an initial condition and w{x,t) is a vector of test functions, which are 
differentiable with respect to x and t as many times as is required and are identically 
zero outside D{x,t). This equation is derived by multiplying the original equations by 
w and integrating the product by parts. Note that differentiability of weak solutions is 
not required. If the solution is differentiable and satisfies Eq. (|2.ip . then Eq. (|2.5p is also 
satisfied. Hence ordinary solutions are indeed weak solutions. It should be remarked that 
uniqueness of weak solutions for a given initial condition is no longer guaranteed. 

2.2. Solutions of Riemann problem: simple waves and discontinuities 

The Riemann problem is initial value problems for a system of conservation equations, 
in which two constant states are separated initially by a discontinuity, which is located, 
e.g., at a; = 0: 

uix,0) = \ ;">^) (2.6) 

^ ^ \ Ml (a; < 0). ^ ^ 

We seek weak solutions to this problem. As remarked at the end of the previous section, 
the weak solutions are not unique in general and it is one of the important issues in 
the Riemann problem to determine which solution is physically meaningful. Since there 
is no typical time and length scales in the problem, the solution should be self-similar. 
It is well known that solutions of the Riemann problem consist of a couple of centered 
simple waves and discon tinuities. The latter satisfies the Rankine-Hugoniot relations (e.g. 



Jeffrey fc Taniutilll964D 



The simple waves are defined as waves, in which there is only one independent com- 
ponent in u and all the other components are its functions. They are related with the 
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Figure 1. A schematic picture of the fc-th Riemann invariants and corresponding locus of the 
fc-th simple wave in the case of = 3. There are two fc-th Riemann invariants: Ji and . The 
red surface is defined by = const, and the blue one is given by = const. The intersection 
is a line called the locus of the fc-th simple wave, to which the fc-th right eigenvector, rk{u), is 
tangential. 

eigenvalues of the Jacobian matrix (j2.2l) in such a way that all quantities in the simple 
wave are constant along each characteristic given by 



The simple wave that corresponds to eigenvalue Xk is referred to as the fc-th simple wave. 

Associated with the fc-th simple wave are so-called the fc-th Riemann invariants de- 
noted by (/ = 1, 2, . . . , — 1), which are constant across the fc-th simple wave. They 
are defined to be the quantities that satisfy the following relation with the fc-th right 
eigenvector, r^: 



The fc-th Riemann invariants define A^ — 1 hypersurfaces in phase space and their inter- 
section is a line, to which the fc-th right eigenvector, r^, is tangential, and is called the 
locus of the fc-th simple wave. See a schematic picture drawn in Fig. [TJ 

The Rankine-Hugoniot relations, which are satisfied by the quantities on both sides of 
a discontinuity, are expressed as 



where s is the velocity of the discontinuity and |A"] :~ Xi — Xq denotes a jump in 
quantity X across the discontinuity and and Xi are the values ahead of and behind 
the discontinuity, respectively. Since this is a system of A^ equations, we obtain a one- 
parameter family of solutions for a state given on one side of the discontinuity. Not all 
solutions are physical as mentioned repeatedly. Rarefaction shock waves for hydrody- 
namical equations are well known examples. Since entropy decreases across such shock 
waves, they are certainly unphysical and rejected. This entropy condition is not sufficient 



X = Xkt + constant. 



(2.7) 



(2.8) 



sH - i/i = 0, 



(2.9) 
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to eliminate all unphysical solutions in general and the so-called evolutionary conditions 
are introduced. 



2.3. Evolutionary conditions 

Intuitively put, the evolutionary conditions require that physical solut i ons should be sta- 
ble against splitting into other waves (see e.g . Falle fc Koniissarovl2001 : Polovin fc Demutskii 
IQQOtlLandau et al. l ll984l : IJeffrev fc Taniutilll964l : lLa>dll957ir To be more precise when a 
small-amplitude wave is injected to the discontinuity, the resultant state, which is deter- 
mined uniquely, is required to still consist of the same type of discontinuity and outgoing 
small-amplitude waves. This is satisfied if the number of characteristics emanating from 
the discontinuity is -I- 1 and if the initial jump, is linearly independent of the 
eigenvectors that correspond to the outgoing characteristics. 



3. Ideal MHD 

In plane symmetry, the ideal MHD equations are given by 

dp d 

+ 



d 



d 



dt 



dx 



d_ 

dt 



dBt 

dt 



d_ 

dx 
d_ 

dx 



de d 
dt dx 



e+p + 



{pVnVt - BnBt) 
{v,iBt - BnVt) 



- BnB ■ V 



0, 
0, 

0, 
0, 

= 0, 



(3.1) 
(3.2) 

(3.3) 
(3.4) 

(3.5) 



where p, jp, v and B are density, pressure, flow velocity and magnetic field, respectively 
([Landau et aLMl984l) . The subscripts n and t indicate the normal component, i.e. x- 
component, and transverse component, i.e. y or z-component, respectively. The total 
energy density is denoted by e = p/(7 — 1) -I- pv^ /2 + B^/2, where the equation of state 
for ideal gas is assumed and 7 is the ratio of specific heats. The normal component of 
magnetic field, i?„, is constant owing to the divergence- free condition. 

3.1. Simple waves 

The eigenvalues of the Jacobian matrix for the system of Eqs. p.ip - p.Sp are 

VnTCf, VnTCA, VnTCs, W„, (3.6) 

where c/, ca and Cg are called the fast, Alfven and slow speeds, respectively, and are 
expressed as 



1 fB' 




(3.7) 



(3.8) 



In the above expressions, a = \fyplp is the acoustic speed. In Eq. p.6|) . the minus (plus) 
sign is applied to the left-going (right-going) waves. The simple waves corresponding to 
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these eigenvalues are referred to as the fast, Alfven, slow and entropy waves, respectively. 
The right eigenvectors for fast and slow waves are given as 



± 



-P 
-IP 



1 - {Cf,slcAY Bn 

Bt 



(3.9) 



in which the new variables, ^/ and ^g, are introduced as follows: 



(3.10) 



These factors are necessary to ensure that the eigenvectors do not vanish for any u. In 
deriving the above expressions of right eigenvectors we assume that u = {p,p, Vn,vt,Bt). 
Since we do not use the right eigenvectors for the Alfven and entropy waves in solving 
the MHD Riemann problem, we do not give their explicit forms here. Note that the 
Alfven and entropy waves are linear whereas the fast and slow waves are neither linear 
nor genuinely non-linear. 

The eigenvalues are degenerate in the following two cases: 



B„ 



0: Cf 




CA = 0, 



B„ 7^ 0, Bt = : Cf = max(a, Ca), Cg = min(a, ca) 



(3.11) 

(3.12) 



In the latter case the limits of the right eigenvectors for fast and slow waves as Bt ^ 
depend on the magnitudes of the acoustic and Alfven speeds. For a > ca we obtain 



(3.13) 



where ej is a unit vector that has the same direction as the transverse magnetic field. 
Note that rj is reduced to the eigenvectors for the rarefaction waves in the ordinary 
hydrodynamics. In the opposite case, i.e. a < ca, we get 



-p 







-IP 







±a 


, a 










±sgn(S„) et 

















-P 









-IP 


a 





«.T — i 
) ' s ^ 


±a 




Tsgn(B„)et 









-^/pet 








(3.14) 



in which rf is reduced to the eigenvectors for the ordinary rarefaction waves in hydro- 
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dynamics. Finally in the case of a = ca, we find 



f 



1 

71 



1 

71 



-IP 

±a 
±aef 



(3.15) 



-P 
-IP 
±a 
Taet 

As mentioned earlier, the right eigenvectors are chosen in our code so that these degen- 
erate cases could be properly handled as the limits of non-degenerate cases. 

In the fast rarefaction wave the magnitude of transverse magnetic field is decreased 
and, as a limiting case, it vanishes behind the so-called switch-off rarefaction wave. Since 
the fast rarefaction wave cannot reverse the direction of the transverse magnetic field, 
the switch-off rarefaction is the end point of the fast rarefaction locus. 

3.2. Discontinuities 

As mentioned earlier, discontinuities are another important element in the solutions of 
Riemann problem. The quantities on both sides of a discontinuity satisfy the Rankine- 
Hugoniot relations, which in ideal MHD are expressed as 



m 



pv 



7-1 





m = const.. 


(3.16) 


Iv] + 




= 0, 


(3.17) 


mlvt] 


-B4Bt 


-0, 


(3.18) 


mlvBtj 


-BnlVt 


= 0, 


(3.19) 


{pm + \ 




)^0, 


(3.20) 



in the rest frame of the discontinuity. In the above expressions, m := poVno = PiVni is 
the mass fiux, v := 1/p is the specific volume and {X) :— {Xq + Xi)/2 stands for the 
arithmetic mean of upstream and downstream quantities. The solutions can be classified 
into the fast, slow or intermediate shock or into the rotational, contact or tangential dis- 
continuity. Details of the classification are found i n the literature (e.g. iJeffrev fc Taniuti 
1964t iPolovin fc Demutskii Il990l: lTorrilhon|[2002i |2003| ). In what follows, we summarize 
only those f eatures that are needed for later discussions. 

Following iTorrilhonI (|2002l l2003f) , we normalize all quantities with those upstream as 



Wo' 



-Bin 



B 



- El 

Po' 


Bt : 


Bn 


Bn 

Vpo' 


Mo 


._ Em 

ao 



(3.21) 
(3.22) 



and employ in the following the dimensionless MHD Rankine-Hugoniot relations, which 
are obtained by substituting Eqs. ((3?2Tt - ((3?22t and ehminating {vt] in Eqs. ((3TT6)) - (13^20)) : 



p- l + jA4^{v- 



A)~BHBt 



Mn 



-jipi - 1) + - m + 1) + \{v - 



-A)=Q 



0. 



(3.23) 
(3.24) 
(3.25) 



Fixing the upstream quantities. A, B and Mq, we solve Eqs. p.23|) - p.25|) and use 
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Eqs. p.2ip - p.22p to obtain v, p and Bt- The other downstream quantities can be calcu- 
lated as 

Vnl = VVnO, (3-26) 

Vt,=Vto±^lBtl (3.27) 

In Eq. p.27p . the plus and minus signs correspond to the left- and right-going disconti- 
nuities, respectively. 

3.2.1. contact, tangential and rotational discontinuities 

The solutions of Eqs. (|3.23p - p.25p that have a vanishing mass flux, i.e. Mq — but 
a non-vanishing normal component of magnetic field, i.e. i? 7^ 0, are called the contact 
discontinuity and satisfy the following relations: 

V — arbitrary, p = 1, Bt = A, (3.28) 

H - 0. (3.29) 

The solutions with Mq = and B — 0, on the other hand, are named the tangential 
discontinuity, for which the following relations hold: 

V = arbitrary, p-l + i(s/ - A^) = 0, (3.30) 

Ivnj = 0, Ivtj = arbitrary. (3.31) 

The solutions with Mq 7^ and B are either a linear wave (£> = 1) or a shock wave 
{v > 1). The former is referred to as the rotational discontinuity, since the transverse 
component of magnetic field rotates, not varying its magnitude during its passage. The 
rotational discontinuities meet the following conditions: 

v = l, p=l, Bt^^A\ Mo=— , (3.32) 

7 

Kl-0, lvt\=±^{Btl (3.33) 

where the plus and minus signs correspond to the left- and right-going waves, respectively. 
The last relation in Eq. p.32p implies the upstream Mach number is equal to the ratio 
of the Alfven velocity to the acoustic speed there. 

All the above discontinuities satisfy the evolutionary conditions except for the rota- 
tional discontinuity in which the transverse magnetic field rotates by 180°. Th e latter is 



sometimes called weakly evolutionary in the literature (jJeffrev fc Taniutilll964[ ). since the 
neighboring rotational discontinuities are all evolutionary. Following this point of view, 
we treat the solutions with weakly evolutionary discontinuities as regular solutions in 
this paper U 

3.2.2. shock waves 

The solutions of Eqs. (|3.23p - (|3.25p . for which Mq ^ and -0 > 1, i.e. matter is com- 
pressed as it passes through the discontinuities, are called shock waves. Their notable fea- 
ture is that magnetic fields are either planar or coplanar. This is apparent from Eq. (|3.24p . 
Indeed, recalling i) > 1, we obtain 

Bt = -^—h —:A. 3.34 



\ This may need reconsideration, however, and will be addressed elsewhere. 
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Substituting Vni — vvna to eliminate v, we also obtain the following relation: 

B, = ^^^A, (3.35) 

which shows immediately that transverse magnetic fields are coplanar if and only if the 
upstream flow velocity is super- Alfvenic whereas the downstream speed is sub-Alfvenic. 

The shocks with planar transverse magnetic fields are either fast or slow shocks, the 
former of which amplifies the magnitude of transverse magnetic fields whereas the lat- 
ter reduces it. The shocks that change the direction of transverse magnetic fields are 
referred to as intermediate shocks. The terminology employed in the literature is rather 
confusing. It is remarked that the in termediate shock i n this paper is different from the 
intermediate discontinuity defined in lJeffrev fc Taniuti In fact, the latter is used 



for the discontinuities that satisfy the evolutionary conditions whereas the intermediate 
shocks in this paper are non-evolutionary as will be evident shortly. 

Recalling ^ ^ c^, we assign 1 to the states with super- fast velocities, 2 to those 
with sub-fast and super-Alfvenic velocities, 3 to those with sub-Alfvenic and super-slow 
velocities and 4 to those with sub-slow velocities in the shock-rest frame. With this 
allocation, the fast shock is denoted by 1 — >■ 2 shock, since the upstream velocity is 
super-fast (state 1) whereas the downstream speed is sub-fast and super-Alfvenic(state 
2). Similarly the slow shock is designated as 3 — 4 shocks. The intermediate shocks 
normally belong to one of the following four types: 1— >-3, 1— ?'4, 2— 7^3 and 2 — 4 
shocks. The 1 — > 3 and 2 — > 4 intermediate shocks are called over-compressive shocks 
and 9 out of 14(= 7x2) characteristics run into these shock waves. To the 1 — 3 
and 2 — >■ 4 shocks converge the fast and Alfven characteristics and the Alfven and slow 
characteristics, respectively. The 1 — >■ 4 intermediate shock is doubly over-compressive 
and 10 characteristics of all types go into the shock. In the case of 2 — 3 shock, only the 
Alfvenic characteristic converges to the shock wave and the number of in- and out-going 
waves are right. However, the other evolutionary condition on the linear independence of 
eigenfunctions is violated and the shock is hence classified as non-evolutionary. In some 
cases, the fiow velocity coincides with one of characteristic velocities. We employ a pair 
of numbers to specify those states: "1,2 ", "2,3" and "3,4" represent those states whose 
flow speed are equal to the fast, Alfven and slow speeds, respectively. The shock wave 
with the upstream velocity being super- Alfvenic and the downstream speed being equal 
to the slow velocity, for example, is designated as the 2 — 3, 4 shock. Of our special 
concern among these intermediate shocks are the so-called switch-on (1 2, 3) and 
switch-off (2, 3 — >■ 4) shocks, the details of which will be given shortly. 



3.2.3. jast and slow loci 

We regard the shock solutions of Eqs. (|3.23|) - p.25|) as functions of the upstream Mach 
number (Mq), normal (i?„ or B) and transverse {Bt or A) components of magnetic field, 
divide them into two families and look into their loci in some detail. Since magnetic fields 
in shock waves are either planar or coplanar as pointed out earlier, in the following we 
assume without loss of generality that magnetic fields are confined in the (x, t/)-plane 
and treat Bt and A (^ 0) as scalar variables. 

Eliminating p and Bt from Eqs. p.23|) - p.25p we obtain the following cubic equation 
for the specific volume, v: 
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7M0 



7-1 
7+1 



7 + 1 



(3.36) 



= 0. 



In the deriving this equation, we used the assumption oi v ^ \. The family of the fast 
shock is the solutions characterized by the feature that matter is compressed and the 
transverse magnetic field is amplified by the passage of shock wave. It is then found that 
this branch of solutions satisfies the inequality Wmin < "O < 1, where the minimum is given 

by 



B' 



t).in = max(^-^,I^j. (3.37) 

The loci of the solutions are shown in Fig. [2] as a function of the upstream Mach number 
(Mo) for some combinations of the upstream normal {B) and transverse {A) components 
of magnetic field. 

As seen in the right panel of Fig. [2l the fast locus is divided into two branches for the 
special case of ^ = 0, i.e. the vanishing upstream transverse magnetic field, with Mach 
numbers satisfying the following inequalities: 



1 



c/o < Mo < 



h + i B^ 

7- 1 7 



7- 



1' 



(3.38) 



where Cfo := Cfo/ao is the normalized fast velocity. One of the branches that generates 
non-vanishing transverse magnetic fields by the shock passage is called the switch-on 
shock branch and the other is referred to as the Euler shock branch, in which the trans- 
verse component of magnetic field remains zero. The post-shock specific volume and 
transverse magnetic field are given by 



2+(7-l)Afg 
(7 + 1)^2 



_B2 
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7- 
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(3.39) 
(3.40) 



respectively. In the above expressions, the first options correspond to the Euler shock 
and the second ones to the switch-on shock. The requirement that the quantity in the 
square root be non-negative gives the inequality p.38p . It is noted that the flow speed 
behind the switch-on shock is equal to the Alfven speed. The switch-on shock is hence 
designated as the 1 — >■ 2, 3 shock and is non-regular. It is also noteworthy that the 
Euler shock is essentially a hydrodynamical shock wave and its locus is extended to 
< c/o, where it is smoothly connected to the slow-shock counterpart. The Euler 
shock is evolutionary except the range within which the switch-on shock branch appears, 
i.e. the range satisfying the inequality p.38p . 

The slow family is characterized by the feature that matter is compressed but the 
transverse magnetic field is reduced and in some cases reversed by the shock passage. 
The slow loci are shown in Fig. [3] as a function of the upstream Mach number (Mq) for a 
number of combinations of the upstream normal (B) and transverse (A) components of 
magnetic field. It is evident that some loci are two- valued as a function of Mo. There is 
a maximum Mach number in each locus, at which the downstream flow speed becomes 
equal to the the slow speed. In some cases the minimum values appear in the transverse 
magnetic fields for given A and B. 
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Figure 2. The fast loci for different combinations of the normal (B) and transverse (^4) com- 
ponent of magnetic field. The left panels: B = 1 and A = 2 (red), 1.5 (green), 1 (blue), 0.5 
(purple) and (light blue). The right panels: 5 = 3 and A — 1 (red), 0.5 (green), 0.05 (blue), 
(purple, switch-on shock) and (light blue, Euler shock). The switch-on shock does not exist 
for B = 1. See the text for details. 



The intermediate shocks are those that give negative downstream transverse magnetic 
fields. The shocks that nullify the transverse magnetic field are called the switch-ofT 
shock (2,3 — ^ 4 shock). It is apparent that the switch-ofF shocks are located at the 
boundary between the regular slow shocks and the intermediate shocks, the fact which 
will be important later. Each locus is terminated at the point that corresponds to a 
rotational discontinuity, which rotates magnetic field by tt radian, i.e. Bt = —A. In the 
limit of i? — 7- or A — > 0, the slow branch vanishes, i.e. there is no solution that satisfies 
Eq. (I3.36P and the inequalities: Cso < Mq < c/o and < -C < 1, except when A = and 
the upstream Alfven speed is larger than the acoustic speed, i.e. ca(:= c^o/oo) > !> in 
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Figure 3. The slow loci for different combinations of the normal (B) and transverse (A) com- 
ponent of magnetic field. The left panels: 5=1 and A — 3 (red), 2 (green) and 1 (blue). The 
vertical black dashes indicate the points, at which Mo — cao and switch-off shocks (2, 3 — >■ 4 
shocks) occur. Since cao is independent of A, the Mach numbers, Mq's, at the points for all loci 
coincide with one another. These points mark the boundary between the regular slow shocks and 
non-regular intermediate shocks. The characters, R, and A, attached to each locus stand for 
the regular slow (3 — >■ 4), 2 — > 4 intermediate and 2 — >■ 3 intermediate shocks, respectively. The 
horizontal dash on each locus shows the point, at which the Mach number reaches its maximum 
on the locus and a 2 — >■ 3, 4 shock occurs. This is the boundary between the 2 — >■ 4 (C"*") and 
2 — >■ 3 (A) intermediate shocks. The right panels: B = 3 and A = 3 (red), 2 (green), 1 (blue) and 
(purple, Euler shocks). The vertical black dashes again give the boundary between the regular 
and non-regular shocks, at which switch-off shocks (2, 3 — >■ 4 shocks) occur. The characters, R, 
C"*" and A, have the same meaning as in the left panels whereas C^"*" and C*"*", which emerge only 
for small A's, represent the 1 — >■ 4 and 1 — >■ 3 intermediate shocks, respectively. The horizontal 
dash on each locus marks again the point, at which the maximum Mach number is reached. On 
the other hand, the two vertical blue dashes on each blue locus indicate the points, at which 
Mo = c/o- A 1,2 — >■ 4 shock occurs at the point closer to the vertical black dashes whereas a 
1,2—^3 shock emerges at the other point. A 1 — 3, 4 shock occurs at the point indicated by 
the horizontal blue dash. Note in passing that the locus vanishes at B = 0. 
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which case the Euler shock branch takes its place. This branch is extended to the regime 
of Mq > c/o and connected smoothly to the fast-shock counterpart as mentioned earlier. 



4. Solving the MHD Riemann problem 

The ba sic for find i ng so lutions of the MHD Riemann problem is essentially the same 
as that of lTorrilhonI (|2002l ). However, we have extended it substantially so that we could 
obtain all types of solutions for arbitrary initial conditions, including non-regular ones. In 
the followi ng we briefly describe the ba sic ideas and the technical details will be published 
elsewhere (jTakahashi fc Yamadall2012r) . 

First we consider the regular solutions that were treated bv iTorrilhon ( 2002 ). It is as- 
sumed that both normal and transverse components of magnetic field are non- vanishing. 
This implies that the degenerate case mentioned in § 13.11 is ignored. As a definition of 
the regular solution, no intermediate shock is involved. It is further assumed that there 
is no switch-off rarefaction wave, either. Since the switch-off rarefaction wave is regular, 
this is another genuine restriction on the solutions that can be handled. Thanks to these 
assumptions, which types of waves are generated in what order are known a priori: the 
fast, Alfven, and slow waves are fanning out in this order on both sides of the contact 
discontinuity. Since each regular wave forms a one-parameter family, the regular solu- 
tions considered now are specified by the seven parameters that are associated with the 
waves involved in the solution. At the contact discontinuity, the six quantities other than 
density are continuous. Hence the six parameters out of seven that correspond to the 
regular waves other than contact discontinuity are determined by the requirement that 
the six quantities other than density be continuous at the contact discontinuity. Then 
the last param eter that specifies the contact discontinuity is automatically obtained. As 
pointed out bv lTorrilhoJ(|2002l ). further simplification arises from the observation that 
only the rotational discontinuities rotate transverse magnetic fields. In fact, the sum of 
rotation angles at the left- and right-going rotational discontinuities is equal to the angle 
that are made by the transverse components of the magnetic fields on both sides of the 
initial discontinuity. Hence we have only to solve one of the rotation angles. 

We thus arrive at the reduced problem, in which the remaining five parameters should 
be determined. The system of equations is expressed symbolically as 

FlirL;<ir; uii^i;uL))) - FUrR, ^'^(V'''; un))) = O, (4.1) 

where ip^, ^/'£, ip^, ipfi^ and ip^ are the five parameters that specify the jumps at the left- 
going fast and slow waves, right-going fast and slow waves, and rotational discontinuities, 
respectively. The initial left and right states are given by ul.r — *^{pl,r,Pl.,r,vl,r, Bil.r)- 

The downstream states of the left- and right-going fast waves are expressed as m{ ^, which 

f 

are functions of the upstream states, Ui^ r, and the parameters, V'^ r- Similarly, tt^ ^ 
represent the downstream states of the left- and right-going rotational discontinuities. 
The five downstream quantities, p, v and \B^\, of the left- and right-going slow waves, 
which should be continuous at the contact discontinuity, are denoted by F]^ j^. These 
equations are numerically solved by the Newton-Raphson method in our code. 

As mentioned already, if the switch-off rarefaction is involved in the solution, the above 
method cann ot be apphed as i t is, s ince we do not know a priori which waves emerge in 
the solution. iTorrilhonI (|2002l . |2003|) neglected such possibilities by assuming that both 
normal and transverse compo nents of magnet i c field s are non- vanishing. This is also the 
case for the cod e develope d bv'And reev et al. ( 2008 ) and available online, which appears 
to be based on TorrilhonI (12002, . ,20031 ). Hence their codes cannot handle all the regular 



16 



K. Takahashi and S. Yamada 



solutions in addition to the non-regular ones, which they did not take into account from 
the beginning. 

On the contrary, our code can treat all these cases: vanishing normal or transverse 
component of magnetic field with a possible existence of the switch-off rarefaction waves 
as well as the emergence of intermediate shocks including the switch-on and switch-off 
shocks. All the intermediate shocks other than switch-on shocks are parametrized as 
slow shocks because they are parts of the slow loci. The switch-on shocks and switch-off 
rarefactions, on the other hand, are treated as a part of the fast shock and rarefaction, 
respectively. As mentioned already, since we do not know a priori which waves are gen- 
erated, we need to try all possible combinations. For example, if the transverse magnetic 
fields are non- vanishing initially, we first attempt to find solutions that do not include a 
fast wave on one side of the contact discontinuity, i.e. the solution with cither a " 1 — > 3" , 
" 1 ^ 4" , " 1 ^ 3, 4" , " 1, 2 ^ 3" , " 1, 2 ^ 4" or " 1, 2 ^ 3, 4" intermediate shock being in- 
cluded. We then try to obtain a solution that does include a fast wave but not a rotational 
discontinuity, i.e. the solution with either a "2 — 3", "2 — > 4" or "2 — 3,4" intermediate 
shock or a switch-off shock. Finally, regular solutions are sought. The tech nical details 
of these procedures will be described in a separate paper (Takahashi fc Yamada. 201Z ). 



5. Results 

In this section, we apply our new code t o two specific init ial conditions and their 



neighbors. One of them is first picked up by iBrio fc Wul (|1988l ) and is known to have a 
non-regular solution with a compound wave. This is one of the initial conditions in the 
MHD Riemann p roblem that are best k nown to developers of numerical MHD codes. As 



demonstrated bv lAndreev et al. ( 2008h there is also a regular solution to this problem. 



The controversy in the literature is why the former is almost always obtained in the 
numerical solutions. After confirming the above regular and non-regular solutions by our 
code, we show that there are actually uncountably many non-regular solutions to this 
problem. By studying the properties of these solutions as well as exploring the neighboring 
solutions, we obtain new insights into the controversy. 

In the second application, we present the initial condition, for which there seems to 
be no regular solution. In fact the non-regular solution we found is a limit of a sequence 
of regular solutions and include a switch-off shock wave. We investigate the neighboring 
solutions in detail and reveal how the solution is approached by various sequences of 
solutions. The result forces us to rethink the interpretations of the regular waves and 
switch-off shocks. 

5.1. Brio & Wu problem: uncountably many non-regular solutions 

We present here the solutions fo r the initial condition in the MHD Riemann problem 
solved numerically bv lBrio fc W^ (|l988h . The left and right states in the initial condition 



are given by 

(PL, PL, VxL, VyL, VzL, ByL, B^l) = (1, 1, 0, 0, 0, 1, 0), (5.1) 
(Pfi, PR: v^R, VyR, v,R, ByR, B,r) = (0.125, 0.1, 0, 0, 0, -1, 0). (5.2) 

The normal component of magnetic field is En = 0.75 and we take fo r the adiabatic index 
7 = 5/3 instead of 7 = 2, the value adopted by Brio & Wu' ('1988'). This is because the 
Riemann invariants are singular at 7 = 2 (e.g. Takahashi & Yamada 2012) and cannot 
be used for check as usually done. Note, however, there is no problem to find solutions 
by our code even in this case. In fact, we confirmed that the solutions are not different 
qualitatively between 7 = 5/3 and 7 = 2. This is a coplanar problem, for which both 
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Figure 4. The regular solution for the lBrio fc Wul (| 19881 ) problem. The waves in the red and blue 
portions that are separated by the contact discontinuity are left- and right-going, respectively. 
The designations FR, R, SS and C represent the fast rarefaction wave, rotational discontinuity, 
slow shock and contact discontinuity, respectively. 



regular and non-regular solutions are expected. It is also noted that this can be regarded 
as a Riemann problem of the reduced MHD system. 

We begin with the regular solution. In Fig.|3]we show various quantities in the solution 
at t = 0.1. Vz and are identically zero and omitted in the figure. The solution consists 
of a fast rarefaction wave (denoted by FR in the figure), rotational discontinuity (R) 
and slow shock (SS) going leftwards and a fast rarefaction wave, slow shock and contact 
discontinuity running in the right direction in addition to a contact discontinuity (C). 
No right-going rotational discontinuity exist because the left-going one has rotated the 
magnetic field by tt radian already. This is the same solution as the one obtained by 
Andreev et aZI (|2008l ). 

Next we present a non-regular solution, various quantities in which are displayed at i = 
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Figure 5. The non-regular solution for the Brio & Wu problem that contains a compound wave. 
The notation is the same as in Fig. U except for the designations SR and IS, which stand for the 
slow rarefaction wave and intermediate shock, respectively. 



0.1 in Fig.m As is evident, this solution contains a c ompound w ave (denoted by 2 — 3, 4 
IS and SR in the figure) and is exactly the solution Brio fc Wu ([l988il obtained in their 
numerical simulations. The compound wave is a 2 — ?> 3, 4 intermediate shock, to which a 
slow rarefaction wave is attached. Other waves that constitute the solution are a left-going 
fast rarefaction wave and a contact discontinuity, slow shock and fast rarefaction wave 
going rightward. Note that the compound wave is responsible for changing the direction 
of magnetic field. It should be also pointed out that the strengths of the waves other 
than the compound wave are different between the regular and non-regular solutions. It 
is the merit of our code that these non-regular s olutions can be handled unlike the one 
developed and made available for public use bv lAndreev et al. ( 2008 ). which can treat 
only regular solutions. 
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Another advantage of our code is that it attempts to find all solutions, trying all 
combinations of constituent waves including the intermediate shocks. In its application 
to the current problem, we found that other non-regular solutions exist, which involve 
various 2 — >■ 3 intermediate shocks. Some of such solutions are shown at i = 0.1 in Fig.[Sl 
The solution given in the top panel resembles the non-regular solution with the compound 
wave presented above. A closer look (see the inset), however, reveals a difference: a 2 — > 3 
intermediate shock is closely followed by a slow rarefaction wave but they are detached. 
The separation of the two waves are more apparent in the middle panel, in which the 
strengths of the intermediate shock and rarefaction wave are also weaker than in the top 
panel. In the bottom panel, the solution does not include a slow rarefaction wave but 
a slow shock instead. The 2—^3 intermediate shock is even weaker in this case. It is 
interesting that the strengths of other waves than the intermediate shock do not change 
very much. It should be emphasized here that these are solutions to the same initial 
condition. 

As a matter of fact, these solutions form a one-parameter family, with the Brio & Wu 
solution with the compound wave at one end and the regular solution with a rotational 
discontinuity at the other end point of the sequence. The parameter that characterizes the 
solution family is a magnitude of jump in the transverse magnetic field across the 2 — >■ 3 
intermediate shock. Since the sequence of solutions is obtained by gradually changing 
this parameter, we can deduce that there are uncountably many solutions in fact. As 
already mentioned, of these non-regular solutions and the unique regular solution it is 
the Brio & Wu non-regular solution that numerical simulations almost always obtain. 
From the point of view of the evolutionary conditions, the regular solution is the only 
solution that contains evolutionary shocks alone. Both 2 3,4 and 2^-3 shocks 
are non-evolutionary. Hence numerical simulations seem to prefer the non-evolutionary 
solution to the evolutionary solution. This is the original controversy. Now we have added 
another question: why is the Brio & Wu solution is singled out from the uncountably 
many non-regular solutions with non-evolutionary shocks? 

Since the evolutionary conditions are related with the neighboring solutions, we con- 
sider the initial conditions that are close to the original one. Both non-coplanar and 
coplanar perturbations are discussed in turn separately. As a non-coplanar initial config- 
uration we take the following: 

(PL, PL, VxL, VyL, WzL, ByL, B^l) = (1, 1, 0, 0, 0, 1, 0), (5.3) 

(PR, PR, VxR, VyR, VzR, ByR, B^r) = (0.125, 0.1, 0, 0, 0, cos 3.0, sin 3.0). (5.4) 

The transverse magnetic field in the right state is rotated slightly with other quantities 
being unchanged. In accord to the evolutionary condition, we found only a regular so- 
lution, which is shown in Fig. [T] The solution includes left- and right-going rotational 
discontinuities, which are responsible for the rotation of magnetic fields. As is evident 
from the comparison (see Fig. |4]), this solution is indeed close to the regular solution to 
the original problem and can be regarded as the pert urbed soluti on. On the other hand, 
no neighboring non-regular solution exists in this case. (|l9886l fal) observed in his dissi- 



pative MHD simulations, however, a new type of non-coplanar shock-like structures that 
do not satisfy the Rankine-Hugoniot relations. They concluded that this time-dependent 
"intermediate shocks" are the neighboring states to the ordinary intermediate shocks. It 
is noted, however, that the time-dependent non-coplanar structures of_thesejiew " inter- 
mediate shocks" de pend upon the dissipation coefficients as shown in IWul (|l9886l fall and 
his succeeding one (|w]ll990l) and may not be realized in ideal MHD. 



It is of great interest to investigate coplanar perturbations, since the coplanarity is 
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Figure 6. Some of the uncountably many non-regular solutions for the Brio & Wu problem 
that contain a 2 — ^ 3 intermediate shock. The notation is the same as in Fig. [S] The insets are 
the close-ups of indicated regions. 



maintained even in numerical simulations unless it is broken explicitly. If the symmetry 
is retained indeed, we are essentially dealing with the reduced MHD system, in which 
the Alfven characteristics do not exist. It sh ould be noted then that th e evolutionary 
conditions will be modified as pointed out by iFalle fc Komissarovl (|200l[ ). For example, 
the 1^-3 and 2 — > 4 shocks become evolutionary and the rotational discontinuities 
become non-evolutionary. On the other hand, the 2 — > 3 shock remains non-evolutionary 
with no converging characteristic. It is noteworthy that the 2 — >■ 3, 4 shock is evolutionary, 
since the slow characteristics are marginally converging and the linear independence of 
the eigenfunctions of outgoing waves and the initial jumps is recovered. As a consequence 
of these changes the Brio & Wu solution should be regarded as a regular solution in the 
reduced system whereas the other solutions with 2 — > 3 shocks are still non-regular even 
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in this system. iFalle fc Komissarovl (|l997l I2OOII ) claimed that this is the reason why 



numerical simulations almost always obtain the Brio & Wu solution. The problem may 
be a bit more complicated, however, as will be demonstrated shortly. 

In Fig. [8] we present the solutions we found with our code for the following initial 
condition: 

(pL, PL, VxL, VyL, VzL, ByL, B^l) = (1, 1, 0, 0, 0, 1, 0), (5.5) 

(pfl, PR, VxR, VyR, v,R, ByR, B,r) - (0.125, 0.1, 0, 0, 0, -0.95, 0). (5.6) 

This time we reduced the transverse magnetic field in the right state slightly, retaining 
the coplanarity. We then found essentially the same types of solutions as in the original 
problem: (1) the solution with a rotational discontinuit}!^, (2) the solution with the 
compound wave consisting of a 2 —> 3,4 shock and slow rarefaction wave, and (3) the 
solutions with a 2 3 shock. As mentioned, the second solution is regular and the first 
and third are now non-regular in the reduced MHD system. The three types of solutions 
actually form a one-parameter solution family with the first and second solutions being 
the end points. As shown in Fig. |9] this is also the case for the following initial condition: 

{pL, PL, V^L, VyL, VzL, ByL, B^l) = (1, 1, 0, 0, 0, 1, 0), (5.7) 

{pR, PR, v^R, VyR, VzR, ByR, Bzr) = (0.125, 0.1, 0, 0, 0, -1.05, 0), (5.8) 

in which the transverse magnetic field in the right state is slightly amplified. 

These results show clearly that the Brio & Wu solution, or the regular solution with the 
compound wave, has a unique neighboring solution of the same type as expected from the 
evolutionary conditions in the reduced system. In the case of the non-regular solutions 
either with the 2 — >■ 3 shock or with the rotational discontinuity, however, the perturbed 
solutions are not uniquely determined, since there are too many outgoing characteris- 
tics. And the existence of uncountably many non-regular solutions is a manifestation 
of this under-determination. Unlike the over-determination and, as a consequence, the 
non-existence of neighboring solution, the meaning of the non-uniqueness of perturbed 
solutions is rather obscure. The results of numerical simulations indicate, however, the 
evolutionary conditions in the reduced system appear to work in selecting the solution. 

5.2. An initial condition with no regular solution 

In the continuous modifications of the Brio & Wu problem, we ran into an interesting 
initial condition, for which there seems to be no regular solution. If true, this is a counter- 
example to the conventional expectation that a unique regular solution always exists. 
This initial condition is truly a special point unlike the Brio & Wu problem. In fact, the 
regular and non-regular solutions of the Brio & Wu problem is shared by the neighboring 
coplanar initial conditions as demonstrated above. This is not the case for the problem 
considered here. 

The initial condition is given by the following: 

{pL, PL, VxL, VyL, VzL, ByL, Bzl) = (1, 1, 0, 0, 0, 1, 0), (5.9) 

{pR, PR, v^R, VyR, VzR, ByR, Bzr) = (0.125, 0.1, 0, 0, 0, -0.329875, 0),(5.10) 

in addition to Bn = 0.75 and 7 — 5/3. This differs from the Brio & Wu problem only on 
the magnitude of the transverse magnetic field {—ByR = -Bcrit — 0.329875) of the right 



f The rotational discontinuity may be an inappropriate term in the reduced system, in which 
the Alfven characteristic ceases to exist. We use it, however, to make clear the correspondence 
to the full system. 
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Figure 7. A non-coplanar neighboring solutions to the Brio & Wu problem. The initial mag- 
nitude of the transverse magnetic field of the right state is the same as the original value. The 
notation is the same as in Fig. 3] 




Figure 8. Some coplanar neighboring solutions to the Brio & Wu problem. The initial magnitude 
of the transverse magnetic field of the right state is reduced from the original value. The notation 
is the same as in Fig. [S] The insets are the close-ups of indicated regions. 




Figure 9. Some coplanar neighboring solutions to the Brio & Wu problem. The initial magni- 
tude of the transverse magnetic field of the right state is increased from the original value. The 
notation is the same as in Fig. [S] The insets are the close-ups of indicated regions. 



Regular & non-regular solutions of the MHD Riemann problem 25 

stat^. This is again a coplanar configuration. The only solution our code obtained is a 
non-regular solution, which is presented in Fig. [10] at t — 0.1. Note that Vz and Bz are 
identically zero in this solution. As can be seen, the solution consists of a fast rarefaction 
wave and switch-off shock, both of which go leftward, and a switch-off rarefaction wave 
and Euler shock, which proceed rightward, as well as a contact discontinuity. Although 
it may seem that the shocks and rarefaction waves form compound waves on both sides 
of the contact discontinuity, they are actually detached. As a matter of fact, it is theo- 
retically shown that the switch-off waves, irrespective of shock or rarefaction, never form 
compound waves. As mentioned earlier, the switch-off shock is a 2,3 — ?> 4 shock, with 
the Alfven speed being equal to the upstream flow speed. Preceded by the switch-off 
rarefaction, which is actually one of the fast rarefaction waves, the fast speed and Alfven 
speeds are degenerate ahead of the Euler shock, which is itself an evolutionary shock 
not only in the full system but also in the reduced MHD system in this case although 
it is possible to be non-evolutionary in principle. It should be stressed that it is the 
advantage of our code that c an tr eat these degenerate cases. In fact, other codes, e.g. 



developed by lAndreev et ali (|2008f ) , ignore the switch-off waves entirely. It is normally 
impossible to prove the non-existence of particular solutions by numerical computations. 
We think, however, that no regular solution exists indeed to the current problem. In 
order to demonstrate why we think so and elucidate the special place that the current 
initial condition occupies, we study the neighborhood of this problem in detail. 

Presented in Fig. [TT] are the solutions to coplanar neighboring problems. More specifi- 
cally, we vary the strength of the transverse magnetic field of the right state, Byji. The 
critical value i?crit is approached both from above and from below. We are concerned 
here only with the regular solutions. As shown in the figures, in the case of \Byfi\ > Scrit 
the solutions consist of a fast rarefaction wave, rotational discontinuity and slow shock 
going leftward and a fast rarefaction wave and slow shock going to the right in addition 
to a contact discontinuity. As the initial strength of the transverse magnetic field in the 
right state becomes smaller, the left-going rotational discontinuity and slow shock come 
closer to each other. Among the right-going waves, the fast rarefaction wave approaches 
the slow shock. At the critical point the left-going rotational discontinuity and slow shock 
merge into the switch-off shock and the right-going fast rarefaction wave and slow shock 
become the switch-off rarefaction wave and the Euler shock, respectively. 

As the initial strength of the transverse magnetic field becomes even smaller and 
< -Bcrit is satisfied, the right-going rotational discontinuity emerges and is respon- 
sible for the inversion of the transverse magnetic field. The Euler shock is replaced by the 
rotational discontinuity and a slow shock whereas the switch-off rarefaction is modified 
to an ordinary fast rarefaction wave again at first as long as |-Byfl| is not very different 
from i?crit- A fast shock appears instead of the fast rarefaction wave when be- 
comes smaller than a certain value. On the left side of the contact discontinuity, on the 
other hand, the switch-off shock is changed to an ordinary slow shock and no rotational 
discontinuity exists. It is stressed that the unique non-regular solution at the critical 
point is connected continuously to the sequences of the regular solutions obtained for 
\ByR\ > -Bclit and \ByR\ < Scrit, respectively. 

Next we turn to non-coplanar neighboring problems. This time we change the direction 
of the initial transverse magnetic field in the right state, keeping its strength at the 
critical value BtR = -Bcrit- In this case we always found regular solutions alone, which 
are displayed in Fig. [12] As is evident, there appear left- and right-going rotational 



^ The ratio of the specific heats is not set to 2, the value adopted by I Brio fc Wul ([mi), but 
to 5/3 for the same reason as given for the first case. 
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discontinuities when the transverse magnetic fields are misahgned initially. They are 
preceded by a fast rarefaction wave and followed by a slow shock running in each direction. 
As the misalignment gets smaller, or the critical point is approached, the amplitudes of 
the transverse magnetic fields rotated by these rotational discontinuities become smaller. 
And they disappear at the critical point and we recover the solution including the switch- 
ofF shock and rarefaction wave presented above. Again this clearly demonstrates that the 
non-regular solution at the critical point is a limit of the sequence of the regular solutions. 
Put another way, there are only regular solutions in the sufficiently close vicinity of the 
non-regular solution at the critical point. This is the main reason why we think that there 
is no regular solution at the critical point. Note, however, that there might be another 
regular solution that are not included in the neighborhood and eluded our search. 

It is interesting to see how a coplanar neighboring solution is approached by non- 
coplanar solutions. Note that they are both regular solutions. As observed above, the 
left- and right-going rotational discontinuities cooperate to rotate the transverse magnetic 
field by the angle imposed by the initial condition. It is found that it is always the right- 
going rotational discontinuity that rotates the transverse magnetic field by a larger angle 
if < Btfl < -Bcrit a-nd vice versa. As the misalignment becomes smaller and the coplanar 
configuration is approached, one of the rotational discontinuities that gives a smaller 
rotation gets weaker and disappears for the coplanar configuration. For a non- vanishing 
but very small misalignment, the relative amplitude of the two rotational discontinuities 
as a function of BtR changes very rapidly near BtR ~ Bent - In this sense the non-coplanar 
solutions shown in Fig. [12] for BtR — Bent are rather special themselves. 

It is finally a legitimate question to ask what has become of the coplanar non-regular 
solutions including 2 — > 3 intermediate shocks that are observed near the Brio & Wu 
problem in the previous section. This is an interesting issue in its own right indeed. 
As long as |-Bj,_r| > -Bdit is satisfied, we find uncountably many non- regular solutions 
that include a 2 — > 3 intermediate shock. The sequence of solutions is terminated at 
one end with the regular solution including a rotational discontinuity. The other end of 
the sequence, on the other hand, is initially the 2 — > 3,4 intermediate shock. As |-Bi,_r| 
becomes smaller, however, the 2 — > 3, 4 intermediate shock cannot be reached at some 
point and the sequence is terminated with a 2 ^ 3 intermediate shock. At the same time, 
another non-regular solution with a 2 — > 4 intermediate shock appears. This solution 
moves up the slow locus as \ByR\ gets smaller, becomes a 1,2 — > 4 intermediate shock 
eventually at ByR k, —Bs — —0.50303 and disappears. In the other case of \ByR\ < Bent, 
there are again uncountably many non-regular solutions, which consist of a right-going 
2 — > 3 intermediate shock and a slow shock, in addition to the regular solutions that have 
right-going rotational discontinuities. The latter is an end point of the sequence of the 
former. No solution with a compound wave is found in this case and the other end point 
of the sequence is the 2 — > 3 intermediate shock, which approaches the regular solution 
as \ByR\ increases and seems to merge at ByR —Bq — —0.25. Thereafter only the 
regular solutions exist. Note finally that all these solutions disappear at ByR = —Bent 
and only the non-regular solution with the switch-off shock waves exists. 



6. Summary and Discussions 

In order to explore non-unique solutions to the MHD Riemann problem, we have devel- 
oped the numerical code to obtain weak solutions that satisfy the entropy condition, the 
minimum requirem ent that the physica lly relevant solutions have to meet. Unlike other 
preceding ones (e.g. Andreev et aLlfioOsh . our code can treat both regular and non-regular 
solutions and can also handle the initial condition, in which no transverse magnetic field 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 


-0.2 

1.0' 
0.8 
0.6 
0.4 
0.2 




Regular & non-regular solutions of the MHD Riemann problem 

I.Or 



27 






Switch-off Shock 




ES 




Switch-off 




Rarefaction 







-0.4 -0.2 


0.2 0.4 


X 




\ 






/switch-off 




/Rarefaction 


Switch-off Shocl< 1 


-0.4 -0.2 


0.2 0.4 




0' 
-0.1 
-0.2 
-0.3 

-0.5 
-0.6 
-0.7 
-0.8 
-0.9 

1.0' 
0.8 
0.6 
By 0.4 
0.2 


-0.2 



\ 


Switch-off 




Rarefactiof 


Switch-off Shocl< 






-0.4 -0.2 C 


) 0.2 0.4 


X 





Switch-off Shocl< 




Switch-off 
Rarefaction 



-0.4 



-0.2 



0.2 



0.4 



Figure 10. The non- regular solution with both a switch-off shock and switch-off rarefaction 
wave. The notation is the same as in Fig. \S\ except for the designation ES, which stands for 
the Euler shock. The initial condition is the same as for the Brio & Wu problem except for the 
strength of the transverse magnetic field of the right state. See ^5.21 of the main body for the 
precise value. 



exists. We have confirmed that the shock waves included in solution satisfy the Rankine- 
Hugoniot jump conditions and the rarefaction waves keep the Riemann invariants con - 
stant indeed. Details of the code will be presented elsewhere ( Takahashi fc Yamadal[2012 ) . 
The code will be useful not only for the study of various solutions as demonstrated in 
this paper but also for numerical simulations, in which the code will provide the exact 
solutions to be compared with and may be implemented to evaluate numerical fluxes. In 
this paper we have applied the new code to two initial conditions, which we believe are 
useful for the consideration of which weak solutions are really physical, the problem with 
a long history of controversy. 

The first problem is one of the best known problems for developers of MHD codes: 
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Figure 11. Coplanar neighboring solutions to the solution presented in Fig. 1101 The red line 
represents the transverse magnetic field whereas the green line displays the density. Only the 
strength of the transverse magnetic field of the right state is varied. See t|5.2l of the main body 
for details. 
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Figure 12. Non-coplanar neighboring solutions to the solution presented in Fig. [10] as well as in 
the bottom right panel in this figure. The green line displays the density whereas the red, purple 
and blue lines represent the transverse magnetic fields, each on a diflerent plane, the angle of 
which is shown by an arc. Only the direction of the initial transverse magnetic field of the right 
state is varied in these computations. See i ]5.2l of the main body for details. 



the Brio & Wu problem (|Brio fc Wulll988l ). This particular problem has been frequently 
utilized for the validation of newly developed numerical codes partly because it has a 
coplanar configuration and 2-dimensional codes can be applied. It has been also known 
for a long time that numerical solutions always involve the compound wave, which is a 
2 — >■ 3, 4 intermediate shock, to which a slow rarefaction wave is attached. This is a non- 
regular solutions, in which the compound wave does not satisfy the evolutionary condition 
and its reality has been a subject of debate. There are two compe ting arguments so far. 
Wu a nd his collaborators (IWijll987l : iBrio fc WijlTflii Imlll988fey . il99(i IWu fc Ken^ 



1992f) insist that the intermediate shocks are stabl e and it is the ev olutionary conditions 
to blame. On the other hand, Falle & Komissarov (jFalle fc Komissaro v 1997, 2001) claim 
that the existence of intermediate shocks in numerical solutions are due to the coplanarity, 
which is inherent in the intermediate shocks. Unless the symmetry is broken somehow, 
the solution will retain the symmetry (coplanarity or planarity) and, as a consequence, 
we are dealing with the reduced MHD system. Then the point is that the 2 — 3,4 
intermediate shock that comprises the compound wave in the numerical solution of the 
Brio fc Wu problem becomes regular whereas t he regular solution tha t contains rotational 
discontinuities, which was shown to exist by lAndreev et ali (|2008[ ) and has been also 
confirmed in this paper, is now non-regular in the reduced MHD system. 

It seems that the existence of the uncountably many other solutions that are non- 
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regula r both in the full and reduced systems supports the claim by Falle & Komis- 
sarov (|Falle fc Komissarovl[l997ll200l[ ). since they have never been realized in numerical 
simulations. It is intriguing to point out, however, that these solutions do not satisfy the 
evolutionary conditions not because the number of outgoing characteristics is wrong but 
because the requirement on the linear independence is not met. This implies that there 
exit too many neighboring solutions to these non-regular solutions and we cannot sin- 
gle out a physically relevant one. This under-determinacy is certainly different from the 
over-determinacy that we commonly find for non-evolutionary s hocks and needs further 
inves t igations. It should be a lso mentioned that some authors (lHadalll994l : iMarkovskiil 
IQQsU l none fc Inutsukal 120071 ) insist that eventually all intermediate shocks satisfy the 



evolutionary conditions if dissipations are taken into account. Although this may ex- 
plain the realization of the compound wave in the solution of the Brio & Wu problem, 
it remains to be addressed why other non-regular solutions that are also evolutionary in 
their analyses but are not produced numerically. They may be linearly unstable. Further 
investigations are certainly needed. 

In the second application we have picked up the initial condition, which appears to 
possess no regular solution. This initial condition was found by continuously modifying 
the original Brio & Wu problem in various ways. In fact, the initial magnetic field has a 
coplanar configuration and is obtained from that in the Brio & Wu problem by reducing 
the strength of the transverse magnetic field in the right state. We have obtained two 
sequences of regular solutions as we approach this initial condition from both directions, 
keeping the coplanarity. They are continuously connected at this point by the non-regular 
solution, which includes both the switch-off shock and rarefaction waves. The switch-off 
shock is non-regular and is located at the boundary between the regular slow shocks 
and intermediate shocks on the slow locus. It is true that we cannot prove the non- 
existence of regular solutions by these numerical calculations but the fact that all the 
neighboring regular solutions, coplanar or non-coplanar, are terminated at this non- 
regular solution suggests strongly that regular solution is unlikely to exist. If true, the 
conventional expectation that there is always a unique regular solution is wrong. We 
are certainly interested in if this non-regular solution is what is realized indeed for this 
initial condition. This will be a difficult task for numerical simulations, however, since 
the initial condition has to be prepared precisely. 

This initial condition is singular also in the following sense: as the coplanar config- 
urations with > -Bcrit approach this critical initial condition, the rotational dis- 
continuity is always left-going whereas the right-going one emerges in the opposite case, 
\Byfi\ < Bcrit, and at the critical condition these rotational discontinuities disappear and 
the switch-off shock and rarefaction wave occur; if the initial transverse magnetic fields 
are misaligned slightly, then the sequence of the regular solutions exhibits a rapid and 
drastic change of the configuration in the vicinity of the critical initial condition: the 
left- and right-going rotational discontinuities exchange their roles. Such sudden changes 
of configuration may have a ramification for the stability of the regular solution for the 
non-coplanar initial configurations with BtR ~ Bcrit- 

We have also found that the sequence of non-regular solutions behaves in an in- 
triguing manner. When we approach the critical initial condition from the regime with 
\ByR\ > i?crit, we always find uncountably many non-regular solutions including a 2 — 3 
intermediate shock, which is terminated at one end with the regular solution that consists 
of a rotational discontinuity. The other end point is initially the 2 — >■ 3, 4 intermediate 
shock that forms a compound wave. For smaller |i3y_R|, however, the end point becomes 
a 2 — > 3 intermediate shock. At the same time another non-regular solution including 
a 2 — > 4 intermediate shock emerges, moving along the slow Hugoniot loci toward a 
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1,2 — >• 4 intermediate shock, at which this type of solution ceases to exist. When the 
critical initial condition is approached from the opposite direction, \Byii\ < Bcrit, we 
again find uncountably many non-regular solutions including a 2 — >■ 3 intermediate shock 
unless \Byii\ is very close to Bait- In contrast to the opposite case, no 2 — 3,4 interme- 
diate shock is realized. The sequence of the non-regular solutions shrinks as \Byii\ gets 
larger and at some point it reduces to the regular solution with a rotational discontinuity 
and no non-regular solution exists any longer. For the critical initial condition even this 
regular solution disappears and only the non-regular solution with the switch-off shock 
exists. These interesting behaviors certainly warrant further investigations. 

After all, it seems that the stability of various intermediate shocks needs major re- 
analysis. The conventional evolutionary conditions are certainly unsatisfactory. Maybe 
the introduction of dissipations and its zero limits should be scrutinized. In so doing our 
new code will be useful to generate all possible solutions. We believe that the survey 
and investigation of other initial conditions and the corresponding solutions, which man- 
ifest some singularities, will provide us with new insights and eventually a clue to the 
understanding of which solutions are physically relevant. 
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